











PARTIAL DIFFERENTIAL EQUATIONS: 
PART-2 


One-dimensional Wave Equation 





Lecture #14 


Superposition principle 


Before we take a closer look at some of the important PDE 
occurring in applications such as wave equations and heat 
equations, let us recall a theorem in linear equations. 


If 4), Ua, ..., Uy, are solutions of a homogeneous linear partial differential 
equation, then the linear combination 


“4 
u= Ley, Cc) U= Cy + Collg tee + Cy, 
c= | 


where the c;,i = 1, 2,...,, are constants, is also a solution. 


Superposition principle is the condition satisfied by a linear system. 


One-dimensional wave equation 





Consider a perfectly flexible elastic string stretched between two 
points at x=0 and x=/ with uniform tension 7. If the string is 
displaced slightly from its initial position of rest and released, 
with the end points remaining fixed, then the string will vibrate. 
The position of any point P in the string will then depend on its 
distance from one end and on the instant in time. Its displacement 
u at any time t can thus be expressed as u=/(x,t) where x 1s its 
distance from the left-hand end. 


One-dimensional wave equation 


The equation of motion is given by 
d-u lod “Ul 
dx” cat’ 


where : c’ = 








vo [A 


in which T is the tension in the string and o the mass per 
unit length of the string. The displacement of the string 
is regarded as small so that T and p remain constant. 


Now let us deal with the solution of this equation. 


One-dimensional wave equation 


fu 1 Ou 


axt = 2 OE has a solution u(x, f). 


The new equation 


Boundary conditions 
(a) The string is fixed at both ends, ie. at x =O and at x =/ for all 
values of time t. Therefore u (x,t) becomes 
eee 9 for all values of t > 0 
u(l, t) =O 


Initial conditions 
(b) If the initial deflection of P at t = 0 is denoted by f(x), then 


u(x, O) = f(x) 
(c) Let the initial velocity of P be g(x), then 


So now we have listed all the information available from the question. 
Next we turn to solving the equation. 


Solution of the wave equation 


Solution by separation of variable. 


Thus, we assume the solution is: u(x,t)=X(x)T(t) or be 
simplified as u=XT. 


ou, fu, 
Ou O-u 
— = XT" d= el 
at ane 6t2 
Ou 1 Ou 


The wave equation can then be written as 


axz 2 At? 


1 
XT — axl” 
xX” 1 TT 
and this can be transposed into ae ae 


Solution of the wave equation 


Notice that the left-hand side expression involves functions of x 
only and that the right-hand side expression involves functions 
of ¢ only. Therefore, if these two expressions are to be equal for 
all values of the separate variables, then both expressions must 
be equal to a constant. 


Denote this arbitrary constant by k. Then we have: 


7, X"—-kX=0 and T”’~—c°kT=0 


Solution of the wave equation 


Let us consider the first of these two equations for different values of k. 
(1) if k=0, XY =0 1. X=a ©. X=axt+b. 
But X=O at x=0 .. b=0 |. X=ax 
“, a=b=0 
and X=O at x=! .. a=0 


.. & =0 which is not oscillatory as the problem requires it to be. 


(2) If k is positive, letk=p* |. X"—p*X=0. 
The auxiliary equation is therefore m*—p*=0O .. m*=p 


“, X =Ae™* + Be 

But X=O at x=O |. O=A+B |. B=-A 

and X=0 at x=1 ». O=Ae*-—Ae -. O= Ale —  P) 

. A=0 .. A=B=0 
Here again X = 0 which is not oscillatory. 
(3) If k is negative, letk=-—p? .. X"+p?X=0. 

The auxiliary/characteristic equation has complex conjugate 
poles at +jp. The solution has the form: 


X = Acos px + Bsin px <n S11) 


which fits the requirements. 


Solution of the wave equation 


The second equation: T'-C*kT =0 
Becomes.... T''+c° pT =0 
And the solution 1s: T =Ccos(cpt) + Dsin(cpt) 


So, the overall solution U=XT 1s: 
u(x,t) = (Acos px + Bsin px) (Ccos cpt + Dsincpt) 
and, if we put cp = » .. p= - this becomes 
u(x,t) = (4 cos x +B sin (C cos At + Dsin At) 


where A, B, C, D are arbitrary constants. 


(2) 


(3) 


Solution of the wave equation 


The result, of course, must satisfy the set of boundary conditions 
which we now turn to. 
(a) uw =O when x = O for all values of t. From this, we get: A=0 
Because, substituting u = 0 and x = 0 in result (3) above 
O=A(CcosrAt+DsindAt) forallt .. A=O 


". u(x,t) =B sin x(C cos At + Dsin At) 


(b) uw =Owhenx=!/ for allt .“. O=B sin “(Cos At + Dsint) 


Now B £ 0 or u(x,t) would be identically zero. ... sin : = 0. 


, Ai = nn where n = | ey as ee A=" for n = 1, 2, 3,... 
Note that we exclude n= 0 since this would also make u(x, t) 


identically zero. 


Solution of the wave equation 


As we can see, there is an infinite set of values of A and each separate 
value gives a particular solution for u(x,t). The values of A are called the 
eigenvalues and each corresponding solution the eigenfunction. 


Putting n = 1, 2, 3,... we therefore have 


nen dx . 
vt) Bain™ (Ces sDsno 


CTT t t 
1 Ai = wa uy = sin 16: cos + Dy sin | 

2C 2 t t 
2 2 == Uy = sin | C2.cos=5™ + Dasin =} 

3 t 
3 A3 = — U3 = sin C: cos—™ + Ds sin 


Solution of the wave equation 


Note that the constant B has been absorbed into the constants C and D 


so that BC = C and BD =e where Cy, Co, C3, ...and Dy, Do, D3, as 
are arbitrary constants. 





. NC 
+ D sin 





nN nczt 
u, = sin "Ice COS 


Since the original wave equation is linear, we have already 
noted that if w=u,, u=u,, u=u;,... are particular solutions, then 
from the superposition principle, a more general solution 1s: 





U=Uu, tu, +U, +... 


Therefore, 








Solution of the wave equation 


We still have to find C, and D, and for this we use the initial 
conditions which we have not yet taken into account. 


(c) Att =O, u(x, O) = f(x) forO<x<!l 


Therefore from (4), u(x, 0) = f(x) = » C,, sin — 


n=l 


(d) Also at t = 0, iF = g(x) forO< x<l 
Ot | 0 


We therefore differentiate (4) with respect to f and put t=O, 
which gives 


2(x)= e, D_nsin 7 
n=1 


Because... 


Solution of the wave equation 








du = NIX ncm . neat NCI ncmt 
=. —— ai) ——¢ 
dt “4 l 


l l ‘  l 
So, with t=O 


“= g(x) = yp, ise ; sin 


Finally we need to involve Fourier series techniques to 
determine C,, and D, so that... 


a3 . na 
a f(x)sin(——x) dx 


5 na { g(x) sin(— x) dx 
NCI * [ 


Solution of one-dimensional wave equation 
SUMMARY 
d-u I d-u 


de? at’ 














| | 2 nx nczt » . NcCazt 
The solution ts: aca ie > sin, = L +2, Sim L } 
n=\ 


ar NIX , a. NIX 
——— x) sin( ——) dx B=—— x) sin(—— dx 
pf fo (—) "Tog |S (—) 


That satisfies the initial and boundary conditions: 


u(O,t) =O 
u(L,t) =O 
u(x,0) = f(x) initial deflection 
- = 2(x) initial velocity 
t 


t=0 


Examples 1-2 
|. Find the solution of the wave equation corresponding to the 


triangular initial deflection f(x): yr 1 
ee O<x<— 

f)=4 4 _ 

And initial velocity g(x)=0. - (L—x) . eae b 


2. Find the solution u(x, t) of the wave equation of the string with 
length L=z , c=1, initial velocity is O, when the initial deflection 1s: 


a. Q.O01sin(3x) . Oe . 
. . Sinusoidal initial deflection 
b. = A[sin(x)-0.5sin(2x)] 


Exercise 
Repeat problem in number 2 with L=1 and c=1, and the initial 
deflection 1s f(x)=kx(/-x). 


l l 
Answer: u(x,t) = BF cos mt SIN 70x + = 3at sin 37x + — Sat sin 57x +...) 


